Orbital-spin ordering in the striped antiferromagnetic state of iron-based 

superconductors 



O 
(N 



O 

o 
I 

Oh- 



O 

o 



in 
o 



X 



Guo-Qiang Liu 

Max- Planck- Institut fiir Festkorperforschung, D-70569 Stuttgart, Germany 
(Dated: January 26, 2013) 

The magnetic properties of iron-arsenides are investigated using the LSDA+f/ approach. In 
addition to one high moment state, we find that a positive U also produces two low moment states 
with m~0.4 fiB and m~1.0 /is. The electronic structures indicate that the low moment states 
originate in the strong orbital hybridization between antiferromagnetic Fe neighbors. Therefore, 
the geometry of the FeAs4 tetrahedron is crucial to the low moment states, which is the key to 
understand the negative pressure coefficient of Tc in LiFeAs. Finally, our theory suggests that the 
superconducting phase is an orbital-spin ordered state, where the orbital spin-moments cancel each 
other. 
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I. INTRODUCTION 

Iron-based superconductors (FeSC) have attracted in- 
tense attention since the discovery of superconductivity 
in LaFeAsOi-^jFa; with Tc of 26 Kii A central topic 
of the studies is to understand the mechanism behind 
superconductivity. Since the electron-phonon couphng 
strength was proved to be insufficient to explain the rel- 
ative high Tc both theoretically and experimentallyrTi^ 
spin fluctuations are expected as an alternative candi- 
date for the pairing glue. Therefore understanding the 
magnetic behavior of the parent compounds is urgently 
needed. 

Unexpectedly, the parent compounds of FeSC exhibit 
a wide variety of magnetic properties despite their sim- 
ilar crystal structures. The undoped compounds are 
antiferromagnetic metals with the moments ordered in 
stripes, i.e. ferromagnetically between the nearest Fe 
neighbors in one direction (x) and antiferromagnetically 
in the other (y). The moment of the 1111 compounds 
RFeAsO (R=rare earth) is ~ 0.3-0.8 ^is'^'- In the 122 
compounds AFe2As2 (A=Ca, Sr, Ba, Eu) it is ~ 0.8-1.0 
yUB)2rii whereas in the 111 family the moment is only 
~ 0.09 fiB in NaFeAs^^ and could not be observed in 
LiFeAsi^'^ In the 11 family, no local moment is found in 
FeScfi^ while FeTe has a unique magnetic structure, a so 
called double-stripe order, with a fairly big moment ~ 
2.0 /ifi.^^ This spread in magnetic moments is difficult to 
understand since the parent compounds have very similar 
lattice structures; they all have similar FeX (X=As, Se, 
Te) layers, and density functional theory (DFT) calcula- 
tions give quite similar paramagnetic band structures for 
all the materials. 

It is conceivable that theoretical studies may clarify 
this puzzle. DFT with the local spin density approxima- 
tion (LSDA) or the generalized gradient approximation 
(GGA) successfully reproduced the single-stripe antifer- 
romagnetic (AFM) order in 1111 and 122 materials J^"— 
However the calculated magnetic moment is about 1.5-2 
fiB, which is much larger than the measured magnetic 
moment. Many theories have been proposed to explain 



the outstanding discrepancy between DFT and experi- 
ment. Yildirim applied a fixed-spin-moment calculation 
to LaFeAsO and found the ground state with a moment 
of TO = 0.48 /i_B.^^ Nakamura et al. applied the LSDA-fU 
approach to LaFeAsOi_j;Fa; and reproduced the exper- 
imental moment with a negative Mazin and Jo- 
hannes argued that the small moment is due to the pres- 
ence of antiphase boundaries.-^i^ Bascones et al. presented 
a five-band model and found a low moment state char- 
acterized by anti-parallel orbital magnetic momentsi^ 
However, it is still not clear which picture is more re- 
alistic. 

In this work, we present an LSDA-I-C/ study to ex- 
plain the puzzle of magnetism in FeSC. We show that 
the measured small magnetic moment can be reproduced 
with positive U, and its origin is a Coulomb-enhanced 
orbital hybridization between antiferromagnetic neigh- 
bors. Since the iron compounds are weak to moderately 
correlated systems, ^'^'^^ the around-mean-field (AMF) 
functional^! is employed in this work. All the calcula- 
tions are performed with the full-potential linear aug- 
mented plane wave (FLAPW) within the local-density 
approximation, as implemented in package WIEN2K.— 
We concentrate on the 1111 material LaFeAsO and 
use the experimental lattice structure^ with the space 
group P4/nmm, a=4. 035 A, c=8.7409 A, ZLa=0.1415 and 
ZAs=0.6415. To study the single-stripe AFM order, we 
use the \/2 x \/2 x 1 supercell containing 4 iron atoms. 
144 k points are used within irreducible Brillouin zone. 
The muffin-tin sphere radii of 2.24, 2.28, 2.02 and 1.98 
ao are used for La, Fe, As and O, respectively. 



II. RESULTS AND DISCUSSION 



A. Magnetism and band structure 
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Since LSDA-I-J7 may give multiple minima;^ we 
change the starting magnetic moment on Fe to search 
for all the solutions. Fig. [T] shows our calculated local 
magnetic moment on Fe and the total energy as a func- 
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FIG. 1: Calculated magnetic moment and total energy for 
LaFeAsO as a function of U. 



tion of the Coulomb interaction U. We used two values 
for the Hund's coupling, J=0 and J=0.82 eV, the latter 
having been obtained from an X-ray experiment. Sur- 
prisingly, two kinds of low moment (LM) in addition to 
one high moment (HM) solutions are found in our calcu- 
lations. The HM state was expected since the LSDA+J7 
is usually believed to have a tendency to enhance the 
magnetic moment. The two LM states have the magnetic 
moments m~1.0 fiB (LMl) and m~0.4 hb (LM2), which 
are in good agreement with the measured moments in 
the 122 compounds and in most of the 1111 compounds, 
respectively. We did the same calculation for the 122 
materials BaFe2As2, and also found two LM states, very 
similar to LaFeAsO. 

As shown in Fig. la and lb, the LMl state only exists 
within a certain U range, and it can jump to the HM or 
LM2 states with a discontinuous change of the magnetic 
moment when U decreases or increases. The LM2 state 
requires larger U and shows slight variation in moment 
when U is increasing. Fig. Ic and Id present the total 
energy of the LMl and LM2 states relative to the HM 
state. The figure shows that a small U favors the HM 
state and a large U favors the LM2 state. Consequently, 
the experimental ground state of LaFeAsO^ can be re- 
produced in our calculation if U is large enough. The 
effect of J may also be seen in Fig. 1. In general, the 
results for J—0 and J=0.82 eV are similar. We therefore 
only consider the case of J=0 in the following parts of 
this work. 

Our theoretical band structures of LaFeAsO for the dif- 
ferent states are presented in Fig. [51 The LMl and LM2 
states show very different band structures. The LMl in- 
dicates a rather good metal while the LM2 state is nearly 
semiconducting. Therefore, the magnetic strength and 
electrical conductivity are correlated properties in the 
iron compounds. The coexistence of the LMl and LM2 
states in our calculations may clarify some seemingly in- 
consistent experiments. McGuire et al^ reported a sys- 
tematic study on the 1111 materials. They found that 
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FIG. 2: Theoretical band structures of LaFeAsO for single- 
stripe AFM state. All the LSDA-|-[/ band structures are for 
J=0 and (7=3.26 eV. 

CeFeAsO has a magnetic moment about 0.3 fis and its 
low temperatures electrical resistivity is po=2 mTl cm. 
While neutron scatting experiment by Zhao et al.^ indi- 
cated a magnetic moment of 0.8 /is for CeFeAsO, and 
Jesche et al^ reported an electrical resistivity of po=140 
/iil cm for CeFeAsO. It should be noticed that McGuire 
et al. used polycrystalline samples and Jesche et al. used 
large single crystals. Moreover, a recent neutron scatting 
experiment^ found single crystal LaFeAsO has a mo- 
ment of 0.8 fiB, which is also very different from the ear- 
lier study.— From our theory, these discrepancies can be 
naturally explained by assuming that the samples were 
in different states, which may be due to the quality of 
samples. 



B. Origin of the low moment states 

In the LSDA+U calculations, the on-site elements of 
the density matrix may contain some important infor- 
mation about the detail of the states. In Table U we 
present the occupancy of the Fe 3d orbitals and the con- 
tribution of each orbital to the spin magnetic moment. 
The main elements of the spin-dependent local density 
matrices can be extracted from this table. The main 
difference between the three LSDA-I-U states lies in the 
moment of the xy orbital. In the HM state, nixy is en- 
hanced from the LSD A value of 0.42 hb to 0.52 hb- But 
in the LMl and LM2 states, rrixy is reduced to and -0.43 
/is! This unusual reduction implies that the xy orbital is 
the key to understand the LM states. The HM, LMl and 
LM2 states may be characterized as mxy>0, rrixy^O and 
mxy<0, respectively. We may therefore name the LM 
states as 'orbital-spin ordered' states. In the LM2 state, 
rrixy and rriyz have opposite signs from m^2_y2 and rrixz- 
It is worth noting that this result is in very good agree- 
ment with the tight-binding study^l by Bascones et al. 

Compared with the LSDA solution, LSDA-I-C/ rear- 
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TABLE I: Occupancy of the Fe 3d orbitals and orbital con- 
tribution to spin magnetic moment. The data for different 
solutions correspond to the band structures in Fig. 2. All the 
LSDA+U data are obtained with J=0 and f/=3.26 eV. 
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ranges the orbital occupancy and orbital spin-moment. 
The LSDA-I-J7 method can be mapped onto a single-band 
picture if there is no strong orbital hybridization. In this 
case, the Coulomb interaction has two possible effects 
on the magnetic moment. The first is to enhance the 
exchange splitting between the majority spin and minor- 
ity spin, which causes larger moment for each orbital. 
The second is to induce a charge transfer from the less 
occupied orbitals to more occupied orbitals, which may 
decrease the moment of the most occupied orbital. As 
may be seen, the HM solution agrees well with this pic- 
ture, but the strong rearrangement of vrixy and 777.3^2 
in the LM solutions cannot be understood. The defi- 
ciency of the single-band picture suggests that interor- 
bital coupling plays an important role in the LM states. 
This has been realized by Bascones et al^ They stressed 
that it is important to include all five d orbitals in model 
studies. 

Fig. [3] shows the band character for the and 
d^xy orbitals, where cr stands for the majority spin and 
for the minority spin. One may see that the hybridiza- 
tion of these two orbitals is strongly enhanced in the LM 
states relative to the LSDA solution. This strong hy- 
bridization between different spin channels (or more pre- 
cisely, between antiferromagnetic neighbors) rationalizes 
the reduction of m^y and 771322 in the LM states. Sim- 
ilar enhanced hybridization is also found between some 
other antiferromagnetic orbitals, such as dZ-, .2-dZ,, and 
dy^-d%y. Here we focus on this d'^^^ _^2-d%y hybridization. 
The strength of orbital hybridization depends on two fac- 
tors: hopping integral and orbital on-site energy. In the 
LSDA-|-f7 approach, U does not change the hopping in- 
tegrals, but shifts the on-site energies. As shown in Fig. 
121 the d'^^2_j.2 and d^y bands are much closer in the LM 
solutions than those in the LSDA solution, which is re- 
sponsible for the enhancement of orbital hybridization. 

In a tight-binding model, the d orbital on-site energy 
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FIG. 3: Band character for the LSDA, LMl and LM2 solu- 
tions. The left panels are for the Zz^ — r^ orbital with majority 
spin, the right panels for the xy orbital with minority spin. 

can be written as 

e^=60±i/777 -C/<-KeSc- (1) 

where e stands for the on-site energy, a for the orbital, 
a for the spin, m for the local magnetic moment, 77 for 
the occupancy number. / is the Stoner exchange inte- 
gral, which is about 0.9 eV in metallic Fcj^ The minus 
(plus) sign of the term ^/777 is for the spin that is par- 
allel (anti-parallel) to the direction of the local magnetic 
moment. eJJ^ corresponds to the double counting term 
of the LSDA-t-J7 method. In this work, we use the AMF 
functionaP^ and have ej^, — U "-q/5 for the d orbitals. 
For two orbitals with opposite spins, we get 

{el-el) = -Im-U{nl~n$) + {el,,-el,) + {^l-el). (2) 

It is reasonable to assume (eg^^-r^ ~ ^Sy) = 0, as the 
LSDA solution shows that the 3z^ — and xy orbitals 
have similar occupancies. Then we can calculate the gap 
between the d'^^2_r2 and d%y orbitals. Using the data 
in Table U we get -1.58 eV for the LSDA, -0.22 eV for 
the LMl and 0.35 eV for the LM2. These values are 
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consistent with the character bands shown in Fig. |3l To 
enhance the hybridization, the d'^^2_j.2 and d'^y orbitals 
need to be close. Eq. [2] shows some competitive factors 
for this requirement. The term —Im opens a gap between 
opposite spins, and therefore the U term and the double 
counting term need to offset this splitting. This gives a 
restriction on U. As shown in Fig. [U there is no LM 
solution when U is too small. For the LSDA state, the 
exchange splitting —Im opens a gap of -1.58 eV. If we 
start the LSDA-I-J7 calculation with the LSDA solution, 
this gap is too large to be overcomed and the system is 
easier to converge to the HM state. To search for the LM 
solutions, we need to start the LSDA+C/ calculation with 
a small magnetic moment. This may be the technical 
reason why the LM states were not found before. 

We have shown that the role of U in the LM states is to 
enhance the orbital hybridization between antiferromag- 
netic neighbors by changing the on-site energies. Here 
we discuss the effect of d-d interorbital hopping, which is 
the other factor influencing the strength of hybridization. 
From the structure of LaFeAsO, we may realize that the 
direct interaction between the d^y and c?322„,.2 orbitals 
is weak, and therefore the hybridization of d'^^2_j.2 and 
d'^y is via the As p orbitals. These p-d interactions, such 
as Pz-dxy and Pz-d^^^-r^, are very sensitive to the Fe- 
As-Fe angle. This characteristic can be used to tune the 
d-d interorbital hoppings. We adjust the Fe-As-Fe angle 
by changing the ratio between lattice constant c and a 
while keeping the volume constant. Our calculated mag- 
netic moment of LaFeAsO as a function of c/a is shown 
in Fig. m The a;-axis of the figure is labelled ZAs/a, 
where zas is the height of the As atom above the nearest 
Fe layer. In the LM2 state, rrixy is anti-parallel to the 
direction of the total moment at small zas/o. as shown 
in TablelH As zas/o- increases, rrixy becomes numerically 
larger and eventually the total moment gets the same di- 
rection as rrixy, which is indicated by the negative sign 
in Fig. m The figure shows that the LSDA state has 
a very weak dependence on the value of zas/o.- Please 
note that the Fe-As distance almost keeps a constant in 
this calculation. On the other hand, the LM states show 
very strong dependence on the FeAs4 geometry. This 
phenomenon confirms the importance of the interorbital 
coupling to the LM states. We may conclude that the 
origin of the LM solutions is the Coulomb-enhanced or- 
bital hybridization between antiferromagnetic neighbors. 
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FIG. 4: Magnetic moment and total energy of LaFeAsO as a 
function of c/a with constant volume. The magnetic moment 
for LMl and LSDA are shifted to fit the window. The experi- 
mental values of zas/o, for LaFeAsO, NaFeAs and LiFeAs are 
indicated by arrows. The LSDA-|-?7 data are calculated with 
J=0 and (7=3.26 eV. 



smaller than the value for LaFeAsO (~2.8 eV). An in- 
teresting question is raised: can LSDA or GGA produce 
the LM states without U7 

Fig. 5 shows the low moment solution we find with 
GGA (without U). The GGA-LM state appears at 
Zyis/a=0.40, and the total magnetic moment reaches a 
maximum at 2;As/a=0-46. The GGA-LM solution is 
not a simple low moment state, but an orbital-spin or- 
dered state. At ZAs/a=0-4:6, the orbital moments are 
TO322_^2=0.02, mj.2_j,2=0.04, mxy=-0.23, mxz=0.53 and 
ruyz— 0.01 p,B- We also present the band structure for 
ZAs/a=0.46 in Fig. 5. The very different band disper- 
sion along X direction and y direction shows that this 
state is far away from the paramagnetic phase although 
it has a small total moment. As may be seen, the GGA- 
LM band structure is a little bit similar to the LM2 band 
structure shown in Fig. 2, while the main difference is the 
Coulomb interaction nearly opens a gap at Fermi level. 
Interestingly, we find that LSDA can hardly produce this 
low moment state. A possible reason is that LSDA con- 
tains less electron correlation in the iron compounds than 
GGA does. 



C. Low moment solution from GGA 

Fig. 1 indicates a critical value Uc for the LSDA-|-t/ 
calculations: the LM states only exist dX U > Uc- 
According to our explanation for the LM states, the 
Coulomb interaction and FeAs4 structure are the two fac- 
tors giving rise to the orbital-spin ordered states. Then 
Uc may be reduced by changing the FeAs4 geometry. We 
find that Uc for NaFeAs is about 1.9 eV, which is much 



The LSDA+C/ method is sometimes questioned due 
to the double counting problem<^i^ The choice of 
LSDA-f-C/ functional is somewhat arbitrary. The GGA- 
LM state proves that the orbital-spin ordering does not 
come from a specific LSDA-|-f7 functional, but originates 
in the tetrahedral structure of FeAs4. The GGA-LM 
state provides a strong evidence for the orbital-spin or- 
dering in FeSC and also confirms our explanation for the 
LM states. 
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FIG. 5: Low moment state in LaFeAsO calculated by GGA 
(WITHOUT U). The left panel is the magnetic moment as a 
function of c/a with constant volume. The right panel is the 
band structure for zas/ci = 0.46, which is corresponding to 
an orbital-spin ordered state with mtot = 0.17, mxy = —0.23 
and rrixz — 0.33 fiB- 



D. FeAs4 geometry, magnetism and high-pressure 
experiments 

It is known that the FeAs4 tetrahedron tends to have 
similar volume in different iron-arsenides, while its ge- 
ometry is tuned by different inserted layers (see Table 
I of Ref. 16). For example, LiFeAs has the largest As 
height zas in all iron- arsenides while its lattice constant 
a is the smallest. Therefore, Fig. |3] allows us to esti- 
mate the magnetic moments of the 111 materials. The 
figure gives a moment of 0.1 for NaFeAs and -0.27 
IJ,B for LiFeAs. Here we consider the LM2 state because 
the LM2 state becomes more and more energetically fa- 
vorable when zas/o- is increasing as shown in Fig. |4l To 
check these estimated values, we did a direct calculation 
using the experimental structures, which gives m=0.04 
fiB for NaFeAs and m=-0.20 fiB for LiFeAs. This agree- 
ment between the estimation and calculation shows that 
the FeAs4 geometry is a very important factor controlling 
the magnetic properties of FeSC, which is responsible for 
the variety of magnetic moments in the 1111 and 111 
materials. Our calculated magnetic moment for NaFeAs, 
TO=0.04 /is, is in good agreement with the measured 
value. But in contrast to the experiment, ^'^ the ground 
state of LiFeAs in our calculation has the single-stripe 
AFM order. Actually, there might be experimental evi- 
dence for the striped AFM order. Borisenko et al. mea- 
sured the Fermi surface of LiFeAs using angle-resolved 
photoemissioni^ They found that the pocket around the 
r point is shaped like a butterfly, which indicates that 
the 4-fold rotation symmetry is broken. This symmetry 
breaking might be due to the striped AFM order. 

Fig. 4 may also explain the behavior of the parent com- 
pounds under pressure. High-pressure experiment shows 
that the undoped LaFeAsO becomes superconducting at 
pressure ~ 2 GPa and reaches a maximum Tc of 21 K 



at ~ 12 GPa. For higher pressures the T^ decreases 
slowly.'^^, LiFeAs shows very different pressure behavior 
from LaFeAsO. LiFeAs is a superconductor without ex- 
ternal pressure. Under pressure, the Tc of LiFeAs de- 
creases linearly at a rate of 1.5 K/GPa,'^'' which means 
that the maximal Tc is at zero pressure. Interestingly, 
a different high-pressure experiment indicated that the 
behavior of NaFeAs is more similar to LaFeAsO than 
LiFeAs,'^* even though NaFeAs and LiFeAs belong to the 
same family. 

To explain these experiments, we need to consider the 
structural distortion under pressure. Pressure experi- 
ments have shown that the Fe-As bond length is very 
rigid while the Fe-Fe distance shrinks much faster when 
pressure increaseSf^MS i.e., zas/cl increased by pressure. 
It is found that chemical doping and pressure can both in- 
duce superconductivity in BaFe2As2 and they have very 
similar influence on the FeAs4 geometry. ^^^'^^ A notice- 
able fact is that chemical doping hardly changes the crys- 
tal volume^i This finding shows that the FeAs4 geome- 
try is the key structural feature to superconductivity;^ 
whereas the volume reduction under pressure is less im- 
portant. Therefore the a;-axis direction of Fig. 3] can be 
regarded as the direction of increasing pressure. 

We consider the LM2 state in Fig. 4. For LaFeAsO and 
NaFeAs, pressure first suppresses the magnetism. After 
the magnetic moment reaches 0, the pressure starts to 
enhance the magnetism. It is known that the emergence 
of superconductivity in the FeSC is related to the sup- 
pression of magnetism. Then the decreasing magnetism 
corresponds to the increasing Tc , conversely, the increas- 
ing magnetism corresponds to the decreasing Tc- And 
the maximal Tc occurs at the magnetic moment m=0. 
For LiFeAs, pressure can only enhance the magnetism as 
shown in Fig. [H Consequently, the Tc of LiFeAs has 
a negative pressure coefficient, dTc/dP < 0. Therefore, 
these different high-pressure experiments can be consis- 
tently explained by our theory. Fig. 2] shows that the 
LM2 state is more and more energetically favorable than 
the LMl state when zas/o, increases; hence the 122 com- 
pounds, which have the LMl ground state at pressure, 
will enter into the LM2 state when pressure increases. 
Then our explanation for the high-pressure experiments 
is also suitable for the 122 materials. Importantly, in this 
picture, the suppression of magnetism is because the or- 
bital spin-moments cancel each other. This suggests the 
superconducting phase is not a 'real' nonmagnetic state 
but an orbital-spin ordered state. 



E. LSDA+(7 functionals 

In this work, we used the AMF LSDA+t/ functional. 
We also checked our calculation with another widely 
used scheme: the fully localized limit (FLL) functional,— 
which is suggested to be more suitable for strongly corre- 
lated systems. We found one LM solution coexists with a 
HM solution in this fimctional. The LM solution is very 
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similar to the LM2 state and shows very similar zas/ol 
dependence. However, we didn't find the LMl solution, 
and FLL energetically favors the HM solution. 

III. SUMMARY 

In summary, we have studied the magnetic properties 
of iron-arsenides using the LSDA+C/ approach. Two low 
moment states are found in our calculation, which repro- 
duce the measured magnetic moments in the 1111 and 
122 materials well. The LM2 state shows good agree- 
ment with the low moment state found by previous tight- 
binding study. Our study reveals that the origin of the 
low magnetization is the strong orbital hybridization be- 
tween antiferromagnetic neighbors. From the picture in- 



dicated above we have learned: Firstly, the geometry of 
FeAs4 is crucial to superconductivity and magnetism be- 
cause it controls the d-d interorbital coupling. Secondly, 
our theory suggests that the superconducting phase is not 
a conventional nonmagnetic state, but a state where the 
orbital spin-moments cancel each other. This implies an 
inherent relationship between the magnetic correlation 
and superconductivity. 
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